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Abstract 

Based on the Riemann- and Caputo definition of the fractional derivative we use the 
fractional extensions of the standard rotation group SO (3) to construct a higher di- 
mensional representation of a fractional rotation group with mixed derivative types. 
An analytic extended symmetric rotor model is derived, which correctly predicts the 
sequence of magic numbers in metal clusters. It is demonstrated, that experimental 
data may be described assuming a sudden change in the fractional derivative pa- 
rameter a which is interpreted as a second order phase transition in the region of 
cluster size with 200 < N < 300. 

Furthermore it is demonstrated, that the four different realizations of higher di- 
mensional fractional rotation groups may successfully be connected to the four fun- 
damental interaction types realized in nature and may be therefore used for a pre- 
diction of magic numbers and binding energies of clusters with gravitational force 
and weak force respectively bound constituents. 

The results presented lead to the conclusion, that mixed fractional derivative 
operators might play a key role for a successful unified theoretical description of all 
four fundamental forces realized in nature. 
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1 Introduction 



Since 1984, an increasing amount of experimental data[T]-[3] confirms an at 
first unexspected shell structure in fermion systems, realized as magic numbers 
in metal clusters. 

The observation of varying binding energy of the valence electron, moving 
freely in a metallic cluster, has initiated the developement of several theoretical 
models. Besides ab initio calculation, the most prominent representatives are 
the jellium modeljl] and, in analogy to methods already in use in nuclear 
physics, phenomenological shell models with modified potential terms like the 
Clemenger-Nilsson model or deformed Woods-Saxon potential[2]-[S] • 

Although these models describe the experimental data with reasonable ac- 
curacy, they do not give a theoretical explanation for the observed sequence 
of magic numbers. Therefore the problem of a theoretical foundation of the 
magic numbers is still an open question. 

A more fundamental understanding of magic numbers found for metal clus- 
ters may be achieved if the underlying corresponding symmetry of the fermion 
many body system is determined. Therefore a group theoretical approach 
seems appropriate. 

In this paper we will determine the symmetry group, which generates the single 
particle spectrum of a metal cluster similar to jellium or phenomenological 
shell models, but includes the magic numbers right from the beginning. 

Our approach is based on group theoretical methods developed within the 
framework of fractional calculus. 

The fractional calculus (Zj-pl)] provides a set of axioms and methods to extend 
the coordinate and corresponding derivative definitions in a reasonable way 
from integer order n to arbitrary order a: 

The concept of fractional calculus has attracted mathematicians since the days 
of Leibniz. In physics, early applications were studies on non local dynamics, 
e.g. anomalous diffusion or fractional Brownian motion [8], [9]. 

In the last decade, remarkable progress has been made in the theory of frac- 
tional wave equations. In 2002, Laskin [11] has derived a fractional Schrodinger 
equation. Raspini [12] has proposed a fractional a = 2/3 Dirac equation and 
has shown, that the corresponding 7 matrices obey an extended Clifford alge- 
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bra, which is directly related to SU(3)-symmetry. We [23] have proposed an 
extension of the ordinary rotation group SO(n) to the fractional case. Based on 
this fractional rotation group, an extended fractional symmetric rotor model 
[21] was presented, which was used for a successful description of the ground 
state band spectra of even-even nuclei. 

The definition of the fractional derivative is not unique, several definitions 
coexist [H]- [22]. Recently [13], the properties of higher dimensional rotation 
groups with mixed Caputo and Riemann type definition of the fractional 
derivative have been investigated. From the four different possible realiza- 
tions of this 9-dimensional fractional rotation group it has been demonstrated, 
that the Caputo-Riemann-Riemann decomposition establishes a fundamental 
dynamic symmetry, which determines the magic numbers for protons and 
neutrons and furthermore describes the ground state properties like binding 
energies and ground state quadrupole deformations of nuclei with reasonable 
accuracy. 

On the basis of this encouraging result, in the present work we will demon- 
trate, that the Caputo- Caputo- Riemann decomposition of the 9-dimensional 
fractional rotation group generates a dynamic symmetry group, which deter- 
mines the magic numbers in metal clusters accurately. Furthermore a com- 
parison with experimental data will lead to the conclusion, that a fractional 
phase transition occures at N > 198. 

Finally we will present arguments, that the four different possible decompo- 
sitions of the 9-dimensional fractional rotation group correspond to the four 
different interaction types realized in nature and will give resaonable predic- 
tions on the sequence of magic numbers in gravitationally and weakly bound 
clusters. 



2 Notation 



We will investigate the spectrum of multi dimensional fractional rotation 
groups for two different definitions of the fractional derivative, namely the 
Riemann- and Caputo fractional derivative. Both types are strongly related. 



Starting with the definition of the fractional Riemann integral 

T a f( ) = ) vr/ " /)(x) = rfe -0^7(0 * > m 



where T(z) denotes the Euler T-function, the fractional Riemann derivative 
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is defined as the result of a fractional integration followed by an ordinary 
differentiation: 



It is explicitly given by: 



nO x J \X) - \ ( x 1 Q rO, CfC \ — a r { £\ l 4 J 



The Caputo definition of a fractional derivative follows an inverted sequence 
of operations ([3]). An ordinary differentiation is followed by a fractional inte- 
gration 



dx 



o% = nl 1 -*^ (5) 



This results in: 



( O a?/)(x) = i ^j?de(x-o- a ^/(0 *>o 



Applied to a function set f(x) = x na using the Riemann fractional derivative 
definition (01) we obtain: 



a^^= r [ 1+nQ) s(-l>° (7) 

r(l + (n — l)a) 

= R [n]x^ a (8) 
where we have introduced the abbreviation R [n]. 

For the Caputo definition of the fractional derivative it follows for the same 
function set: 



q« not 
c u x x 



r(l+(n-l)a) x n ^ u 

n = 

(9) 

[n] x {n ~ 1)a (10) 



where we have introduced the abbreviation n [nl 
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Both derivative definitions only differ in the case n = 0: 



oM=»[b]-<W[0] (11) 

-»W-«-rJlbj (12) 

where 5 mn denotes the Kronecker-5. We will rewrite equations (jSJ) and ffTOl) 
simultaneously, introducing the short hand notation 

^^^cH^ 1 '" (13) 



We now introduce the fractional angular momentum operators or generators of 
infinitesimal rotations in the i,j plane on the A-dimensional Euclidean space: 

R , C L,» = ih(x« R , c d? - x^ofif) (14) 



which result from canonical quantization of the classical angular momentum 
definition (for details see [21] )• The commutation relations of the fractional 
angular momentum operators are isomorphic to the fractional extension of the 
rotational group SO(N) 

r,c [Lij(a>), Lki(a)} = ih K ,cfijki mr ' ''{&) (a) (15) 

k, l,m,n = 1, 2, .., N 

with structure coefficients K,cfijki mn (c()- Their explicit form depends on the 
function set the fractional angular momentum operators act on and on the 
fractional derivative type used. 

The structure coefficients and Casimir-operators of the fractional rotation 
group K SO a (3) based on the Riemann fractional derivative definition have 
been derived in [24J and for c SO a (3) based on the Caputo fractional deriva- 
tive definition are given in [23J. We summarize the major results: 

According to the group chain 

n, c SO a (3) D R , c SO a (2) (16) 

there are two Casimir-operators Aj, namely A2 = L z (a) = Li2(a) and A3 = 
L 2 (a) = L\ 2 (o>) + Lf 3 (a) + L\^{a). We introduce the two quantum numbers 
L and M, which completely determine the eigenfunctions \LM >. It follows 

HiC L z (a)\LM > = ftsign(M) R , c [|M|] \LM > (17) 
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Fig. 1. Spectrum of the Casimir operator L 2 (L,a) from (|18|) as a function of the 
fractional derivative coefficient a. Only the L = state differs for Riemann and 
Caputo derivative. 

M = -L,-L + 1,...,±0, ...,L 

R ,c L 2 {a)\LM > = R , c h\ c [L] H ,, C [L+ 1] \LM > (18) 

L = 0,l,2,... 

where \M\ denotes the absolute value of M. In addition, on the set of eigen- 
functions \LM >, the parity operator IT is diagonal and has the eigenvalues 

U\LM >= (-1) L \LM > (19) 



In figured] the eigenvalues of the Casimir-operator L 2 are shown as a function 
of the fractional derivative coefficient a. As has been pointed out in [21], near 
a ~ 1 there is a region of rotational type of spectrum, while for a ~ 1/2, 
the levels are nearly equidistant, which corresponds to a vibrational type of 
spectrum. 

In addition, for decreasing a < 1 higher angular momenta are lowered. 

Only in the case L = the spectra differ for the Riemann- and Caputo deriva- 
tive. While for the Caputo derivative 

c L 2 (a)|00>=0 (20) 



because c [0] = 0, using the Riemann derivative for a ^ 1 there is a non 
vanishing contribution 



t L a (a) |00 >= fr 2 R [0] R [l]|00 >= h 2 T y +a |00 > (21) 

1 (1 — a) 
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3 The Caputo-Caputo-Riemann symmetric rotor 



We now use group theoretical methods to construct higher dimensional rep- 
resentations of the fractional rotation groups RtC SO a (3). 

For a 9-dimensional rotation group G there exist four different decompositions 
with the following chain of sub algebras: 



R.RR.C D n SO a (3) D K SO a (3) D n SO a (3) (22) 

cr R G D c SO a (3) d R SO a (3) D n SO a (3) (23) 

CCR G D c SO a (3) D c SO Q (3) D R SO a (3) (24) 

CCC G D c 50 a (3) D c SO Q (3) D c SO Q (3) (25) 

As an example of physical relevance in [13] the properties of the 9 dimensional 
fractional rotation group C nnG have already been investigated. The main result 
was the verification of magic numbers for nuclei: 

^nuclear 1 = \( N + V( N + 2 )( N + 3 ) N = 0,1, 2, 3, ... (26) 
= 2,8,20,40,70,112,... (27) 

^nuclear 2 = \( N + + 2)(iV + 3) + 2(N + 2) (28) 
= 6, 14, 28, 50, 82, 126, 184, 258, ... (29) 



In this section we will discuss the properties of the group ccrG'. 

For that purpose, we associate a Hamiltonian H, which can now be written 
in terms of the Casimir operators of the algebras appearing in the chain and 
can be analytically diagonalized in the corresponding basis. The Hamiltonian 
is explicitly given as: 

H=^ c Ll(a) + ^ c Ll(a) + ^ R Ll(a) (30) 

with the free parameters 001,002,^3 and the basis is \L1M1L2M2L3M3 >. Fur- 
thermore, we demand the following symmetries: 

First, the wave functions should be invariant under parity transformations, 
which according to (fl9|) leads to the conditions 

Li — 2m L 2 = 2n 2 L 3 = 2n 3 , n x , n 2 , n 3 = 0, 1, 2, 3, ... (31) 
second, we require 
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a L gl (a)\L 1 M 1 L 2 M 2 L 3 M 3 > = +h c [L l }\L 1 M 1 L 2 M 2 L 3 M 3 > (32) 
c L 22 (a)|L 1 M 1 L 2 M 2 L 3 M 3 > = +h c [L 2 ]\L 1 M 1 L 2 M 2 L 3 M 3 > (33) 
n L Z3 (a)\L 1 MiL 2 M 2 L 3 M 3 > = +K[L 3 )\L X M X L 2 M 2 L 3 M 3 > (34) 

which leads to the conditions 

M 1 = 2m M 2 = 2n 2 M 3 = 2n 3 , m, n 2 , n 3 = 0, 1, 2, 3, ... (35) 

and reduces the multiplicity of a given |2niMi2n 2 M 2 2n 3 M 3 > set to 1. 
With these conditions, the eigenvalues of the Hamiltonian (1301) are given as 

E(a) = Huj 1c [2ni] [2ni + 1] + Huj 2 c [2n 2 ] c [2n 2 + 1] + 

^ 3R [2n 3 ] R [2n 3 + l] (36) 

= ±^ r (l + (2n i + l)a) 
ti t T(l + (2n l -l)a) 

r(l + q) r(l + q) 

-O m0 «WiTvl r-d„ 2 o^2^ r (37) 

1 (1 — a) L (1 — a) 

ni,n 2 ,n 3 = 0, 1,2, .. 

on a basis |2ni2ni2n 2 2n 2 2n 3 2n 3 >. 

This is the major result of our derivation. We call this model the Caputo- 
Caputo-Riemann symmetric rotor. What makes this model remarkable is its 
behaviour in the vibrational region near the semi- derivative a = 1/2. 

On the left of figure [5] we have plotted the energy levels in the vicinity of 
awl/2 for the case 

uji = uj 2 = uj 3 = too (38) 
which we denote as the spherical case. 

For the idealized spherical case a = 1/2, using the relation T(l + z) — zY{z) 
the level spectrum ( 1771) is simply given by: 

3 1 1 

E(a = 1/2) = ftw (ni + n 2 + n 3 + - - -S ni0 - -<5„ 2 o) (39) 

For rii 7^ and n 2 ^ this is the well known spectrum of the 3-dimensional 
harmonic oscillator. Assuming a twofold spin degeneracy of the energy levels, 
we introduce the quantum number iV as 

N = m + n 2 + n 3 (40) 
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Consequently we obtain the multiplets of the 3- dimensional harmonic oscilla- 
tor with magic numbers for 



n HO = -(N+l)(N + 2)(N + 3) N = 0, 1, 2, 3, ... (41) 

= 2, 8, 20, 40, 70, 112, 168, 240, ... (42) 

at energies 

E(N) so = huj (N + 3/2) (43) 



In order to determine the multiplets of ( 17H|) . we distinguish two different sets 
of magic numbers: 

For rii — and n-i = the multiplicity of a harmonic oscillator shell for N 
at energy ( fFTI) is increased by exactly one state, the |0000iV + 1 iV + 1 > 
state, which originates from the N + 1 shell. Therefore we obtain a first set 
ccR^niagic i of magic numbers (including the state 1 000000 >): 

ccR^magic 1 = "ho + 2 N = -1, 0, 1, 2, 3, ... (44) 

= ^(iV+l)(iV + 2)(iV + 3) + 2 (45) 

2, 4, 10, 22, 42, 72, 114, 170, 242, 332, 442, ... (46) 

at energies 

E(N) CCR n magic i = hu Q (N + 3/2) (47) 

In addition, for n x = 0, which corresponds to the 1 00 2n 2 2n2 In^n^ > states 
and tt-2 = respectively, which corresponds to the \ln\lri\ 00 2n^2n^, > states 
with a J2n=i 2 = 2(N + l)-fold multiplicity for each each set rii = and 
n 2 = we obtain a second set ccR^ ma p;i c 2 °f ma gic numbers 



ccR^magic 2 = ccR« mag i c i + 4(iV+ 1) N — 0,1, 2, 3, ... (48) 

= i(iV + l)(iV + 2)(jV + 3) + 2 + 4(iV+l) (49) 

= 8, 18, 34, 58, 92, 138, 198, 274, ... (50) 

at energies 

^(^ccR^maKic 2 = ^o(N + 1) (51) 
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CCR^ 

spherical 



o 

0.44 0.46 0.48 0.50 0.52 0.54 0.56 0.58 0.60 0.62 0.64 0.66 0.1 

a 



Fig. 2. level diagram of the Caputo-Caputo-Riemann symmetric rotor for the spher- 
ical case in the vicinity of a ~ 1/2 

In figure [5], the single particle levels are plotted. A remarkable feature is the 
dominant influence of the |0000iViV > state. For a < 0.5 the harmonic oscilla- 
tor type magic numbers die out. As a consequence, for a ~ 0.48 there remain 
the set of magic numbers ccR n ma gi c 2 shifted by 1, which leads to the series 
2, 4, 9, 19, 35, 59, 93, 139, .... For a « 0.46 the |0000iViV > state has completely 
reached the ccR%n a gic 2 m uftiplet, which in a stable series of magic numbers 
at 2, 4, 10, 20, 36, 60, 94, 140, 200, .... 

On the other hand for a > 0.55 the levels are rearranged to form the set of 
magic numbers of the harmonic oscillator. 

We conclude, that the Caputo-Caputo-Riemann symmetric rotor predicts a 
well defined set of magic numbers. This set is a direct consequence of the un- 
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derlying dynamic symmetries of the three fractional rotation groups involved. 
It is indeed remarkable, that the same set of magic numbers is realized in 
nature as electronic magic numbers in metal clusters. 

In the next section we will demonstrate, that the proposed analytical model 
is an appropriate tool to describe the shell correction contribution to the total 
binding energy of metal clusters. 



4 Binding energy of metal clusters 



We will use the Caputo-Caputo-Riemann symmetric rotor (I77j) as a dynamic 
shell model for a description of the microscopic part of the total energy binding 
energy i?tot of the metal cluster. 



-^tot -^"macroscopic + -^microscopic (^ 2 ) 
= -^macroscopic + ^ (^ 3 ) 

where SU denotes the shell-correction contributions. 



To make our argumentation as clear as possible, we will restrict our inves- 
tigation to the spherical configuration, which will allow to discuss the main 
features of the proposed model in a simple context. We will compare our 
results with calculations for the most prominent metal cluster, the sodium 
(Na) cluster. From experimental data[T],[3], the following sequence of magic 
numbers is deduced: 



"magic Na = 2 ' 8 ' 20 > 40 > 58 > 92 ' 138 > 198 ± 2 > 263 ± 5 ' 341 ± 5 ' 
443 ± 5, 557 ± 5, 700 ± 15, 840 ± 15, 1040 ± 20, 
1220 ± 20, 1430 ± 20 (54) 

For a graphical representation of the experimental magic numbers we intro- 
duce the two quantities: 



6 = iV/(n magic Na (iV + 1) - n magic Na (iV)) (55) 

u = "magic Na( iV + X ) " "magic Na( iV ) ( 56 ) 

Interpreting G as a moment of inertia and u a rotational frequency, figure [3] is 
a backbending plot of the experimental magic numbers. We distinguish three 
different regions of magic numbers. For iV < 200 the plot shows a typical rotor 
spectrum. In the region 200 < N < 300 a typical backbending phenomenon is 
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Fig. 3. Backbending plots of experimentally determined magic numbers for (Na)jv 
clusters from [3], ground state excitation spectrum for }q20s76 an d ground state 
excitation spectrum of {09W74. Squares indicate the experimental values. The se- 
quence of experimentally observed magic numbers may be categorized for N < 200 
to be equivalent to an excitation spectrum of purely rotational type, the thick line 
indicating a fit with a fractional rigid rotor spectrum, for 200 < N < 300 a region 
of backbending type (compare with the plotted ground state band spectrum for 
}5!Os 76 ) and finally, for N > 300 a region with almost constant 9 and a behaviour 
similar to excitation spectra of ug-nuclei (compare with the plotted excitation spec- 
trum of 109 W74). In the upper right the corresponding proposed a(N) from (|6ip is 
plotted. 

observed. For illustrative purposes in figure [3] the same phenomenon is docu- 
mented within the framework of nuclear physics for the ground state rotation 
spectrum of For N > 300 the moment of inertia becomes nearly con- 

stant and the graph may be compared with the rotational K = ^ band of the 
ug- nucleus 109 W74, which is a typical example of a core plus single particle 
motion in nuclear physics. 

These different structures in the sequence of electric magic numbers are re- 
flected in the choice of the fractional derivative coefficient a. For iV < 200, 
a shows a simple behaviour similar to the case of magic nucleon numbers, it 
varies in the vicinity of a ~ 1/2. For the special case of sodium clusters, the 
lowest four magic numbers are reproduced with a > 1/2, while up to iV = 198 
a < 1/2 is sufficient. Within the backbending region there is a sudden change 
in a, which we call a fractional second order phase transition, followed by a 
linear increase of the a value for larger cluster sizes. 

The resulting dependence a(N) is shown in figure [31 
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123/156789 10 11 

N 1/3 

Fig. 4. Calculated shell correction dU for the Caputo-Caputo-Riemann fractional 
symmetric spherical rotor with parameter set (|57p -(|63p 

In order to compare our calculated shell correction with published results, we 
use the Strutinsky method [25j,[26j with the following parameters: 



hw = 


3.96N 


U)x = 


1 


UJ 2 = 


1 


u 3 = 


1 



' [keV] 



a 



7 : 

m ■ 



0.55 N < 43 

0.908 - 0.000834 N N < 51 

0.482 + 0.000025 N N < 260 
0.069 + 0.000139 (N - 260) < 900 

I 0.062 + 0.000083 (N - 260) N > 900 

Ll^oU^ + lJ-HtiV 1 / 3 ]) 3 
4 



(57) 
(58) 
(59) 
(60) 



(61) 



(62) 
(63) 



(1621) follows from the plateau condition dU /dj = and (1631 is the order of 
included Hermite polynomials for the Strutinsky shell correction method. 

In figure H] the resulting shell correction SU is plotted. Magic numbers are 
reproduced correctly within the experimental errors. Furthermore we obtain 
a nearly quantitative agreement with published results for the shell correction 
term obtained e.g. with the spherical Woods-Saxon potential[2]- 



Summarizing the results presented so far, the proposed Caputo-Caputo-Riemann 
symmetric rotor describes the magic numbers and microscopic part of the total 
binding energy for metal clusters with reasonable accuracy. We have demon- 
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strated, that the cluster shell structure may indeed be successfully described 
on the basis of a dynamical symmetry model. 



5 Magic numbers for clusters bound by weak and gravitational 
forces 



We have demonstrated, that the Caputo-Caputo-Riemann rotor correctly de- 
termines the magic numbers of metal clusters and that the Caputo-Riemann- 
Riemann rotor is an appropriate tool to describe the ground state properties 
of nuclei with reasonable accuracy. Both models differ only in the mixing ratio 
of fractional derivatives. The phenomena described, differ only in the inter- 
action type of the constituents which build up the cluster. The behaviour of 
metallic clusters is dominated by electro-magnetic forces, while in nuclei the 
long range part of strong forces is important. 

Therefore we postulate, that the group decomposition 

RRR G d R SO a (3) D R SO Q (3) D R SO a (3) (64) 



will determine the magic number of a cluster, which is dominated by a gravi- 
tational type of interaction between its constituents. 

The Hamiltonian RRR H 

hkk# = p^W + f nLl(a) + ^Lj(a) (65) 



with the free deformation parameters ui, a> 2 , ^3 on a basis IL1M1L2M2L3M3 > 
may be diagonalized and with the symmetries fl3Tj) and fl35|) a levelspectrum 



E(a) = hw x r [271i] k [2tii + 1] + hu 2 R [2n 2 ] R [2n 2 + 1] + 

^3 R [2n 3 ] R [2n 3 + l] (66) 
3 r(l + (2n, + l)a) 

=S^r(i + (^-i)a) (67) 

ni,n 2 ,n 3 = 0, 1,2, .. 
on a basis |2ni2ni2ri22n22n32n3 > results. 

For the idealized spherical case a = 1/2 this spectrum is simply given by: 
E(a = 1/2) = hj ( ni + n 2 + n 3 + ^) (68) 
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which is the spectrum of the deformed harmonic oscillator. In the spherical 
case, magic numbers are determined by: 



rw = ~(iV + 1)(JV + 2)(N + 3) N = 0, 1, 2, 3, ... (69) 
= 2, 8, 20, 40, 70, 112, 168, 240, ... (70) 

at energies 

E(N) RRR = fku (N + 3/2) (71) 

This result may be compared with solutions for an independent particle shell 
model, where the potential is determined by a uniformly distributed gravita- 
tional charge (mass) distribution p(r) = q/V inside a sphere. This potential 
is given by 



v ^ -III F=7i dV < 72 » 



r<R « (73) 



and leads to a radial Schrodinger equation for the harmonic oscillator. 

Therefore we are led to the conclustion, that for microscopic clusters with 
gravitational type of interaction of the constituents there will be variations in 
the binding energy per mass unit according to (1691 . 



Consequently we are left with the fourth decomposition of the 9-dimensional 
fractional rotation group 

CCC G D c SO a (3) D c SO a (3) D c SO a (3) (74) 



The Hamiltonian ccc H 

cccH = ^oL 2 M) + | C L 2 » + ^ c Ll(a) (75) 



with the free deformation parameters Ui, u>2, 003 on a basis \LiMiL 2 M 2 L 3 M 3 > 
may be diagonalized and with the symmetries fl3T|) and fl35|) a levelspectrum 



E{a) = fku lc \2nx\o [2n x + 1] + fku 2 c [2n 2 ] c [2n 2 + 1] + 

hu 3c [2n 3 ] c [2n 3 + l] (76) 



15 



ccc^ 

spherical 
2 2 2 

_J i I i I i I i I i I i I I L_ 



0.44 0.46 0.48 0.50 0.52 0.54 0.56 0.58 0.60 0.62 0.64 0.66 0.68 

a 

Fig. 5. level diagram of the Caputo-Caputo-Caputo symmetric rotor for the spherical 
case in the vicinity of a ~ 1/2 



s r(i + (2m + i)a) 



-5 ni0 hxji 



r(l + a) 

r(i-a) 



r(l + a) 



r(i-a 

ni,n 2) n 3 



T(l + a) 
T(l-a) 



(77) 



0,1,2, 



on a basis |2n 1 2ni2n 2 2n22n 3 2n3 >. results. We call this model the Caputo- 
Caputo-Caputo symmetric rotor. For the idealized spherical case a = 1/2 this 
spectrum is simply given by: 



3 1 1 1 

E(a = 1/2) = ftw (ni + n 2 + n 3 + - - -5 ni0 - -5„ 2 o - -5 n3 o) (7J 
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We obtain a first set ccc% a m C \ of magic numbers 



ccc^magic 1 ~~ n Eo + 6 N — 0, 1, 2, 3, ... (79) 

= i(iV+l)(iV + 2)(iV + 3) + 6 (80) 

8, 14, 26, 46, 76, 118, 174, 246, 336, 446, ... (81) 

at energies 

E(N) ccc n magic i = ^ (iV + 3/2) (82) 
In addition we obtain a second set ccc^ m agic 2 °f ma gic numbers 

ccc^rnagic 2 = ccc^magic 1 = 1? 2, 3, ... (83) 

= i(iV+l)(7V + 2)(iV + 3) + 6(iV + l) (84) 

= 20, 38, 64, 100, 148, 210, 288, ... (85) 

at energies 

E(N) ccc n ic 2 = hu (N + 1) (86) 



Finally, the state 1 000000 > with a 2-fold multiplicity has energy E = and 
therefore does fit into one of the two derived series. 

Consequently we are led to the conclustion, that for microscopic clusters with 
weak interaction type of the constituents there will be variations in the binding 
energy per charge unit according to (1831) . 

Summarizing the results presented in this section, we have associated the 
four different decompositions fl22|) - fl25l) of the 9-dimensional mixed fractional 
rotation group with the four fundamental types of interaction found in nature. 
We found common aspects determining the magic numbers for each group. 
There are always two different sets of magic numbers, one set is a shifted 
harmonic oscillator set, the other set is specific to the group considered. For 
the spherical and idealized case a = 1/2 the four different sequences of magic 
number sets are simply the result of the presence or absence of a Kronecker- 
delta. 

Our investigations lead to the conclusion, that mixed fractional derivative type 
field theories may play an important role in a unified theory including all four 
fundamental interactions. 
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6 Conclusion 



Based on the Riemann- and Caputo definition of the fractional derivative 
we used the fractional extensions of the standard rotation group 5*0(3) to 
construct a higher dimensional representation of a fractional rotation group 
with mixed derivative types, model, which predicts the sequence of electronic 
magic numbers in metal clusters accurately. In the region 200 < N < 300 
we deduced a sudden change in the fractional derivative coefficient a, which 
we interpreted as a second order fractional phase transition. Furthermore we 
have shown, that the microscopic part of the binding energy can be reproduced 
correctly within the framework of this model. 

Hence we have demonstrated, that a dynamic symmetry, generated by mixed 
fractional type rotation groups is indeed realized in nature for electro-magnetic 
and nuclear forces. Within this framework it was possible to make predictions 
for magic numbers in gravity and weak force dominated clusters. 
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